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Abstract 

A simplified model is proposed based on the characteristics of the photomultiplier 
tube (PMT). The Model is compared to the available data, and it successes to 
reproduce the saturation profile of the PMT. The model also adds clarification 
about how the PMT works. 
PACS: 85.60.Ha 

Key words: Photomultiplier; PMT; Double Chooz; R1408; R5912; MCP-PMT; 
PMT Model; Gain; Anode; Dynode 



1 Introduction 



At Drexel University, I was testing the 8-inch, 13 stages, R1408 PMTs which 
will be used for the inner vetos of the Double Chooz experiment [1]. These 
PMTs were used in the 1MB [2] and SuperKamiokande [3] experiments. This 
testing allowed a better understanding of their working. In this paper I propose 
a model inspired from observations during the testing period. The first aim 
of this model is to explain in simple way the saturation characteristics of the 
PMT. The Model gives a global description of how the different PMTs charac- 
teristics, like the quantum efficiency (QE) or the collectivity efficiency (CE), 
are linked together to the PMT gain (G). The model is compared to the data 
obtained from three different PMTs, which have different single photoelectron 
(spe) distributions. 



2 The general case 



The basic idea is at a fixed applied high voltage, the PMT works as charge 
amplifier, which will amplify below a certain charge limit qsat and beyond this 
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limit the PMT will continue to amplify the signal but at lower gain. As the 
signal in the input of the PMT is increased, the PMT gain decreases. 



A general profile which matches this description of the gain variation is given 
by: 

^ , (1) 



Where a is a real coefficient which can be determined from the comparison 
with the data, q is the collected charge at the anode and Gq is the gain at the 
limit as — > , and depends only on the applied high voltage like following 
[4]: 

Go = kV'^''. (2) 



Where k is a constant, which depends on the material of the dynodes and 
the voltage division between them, N is the number of the dynodes and /? is 
between 0.7 and 0.8. Gq is usually set up by appropriate high voltage at low 
light level (single photoelectron level). The Qsat is an intrinsic constant to each 
PMT, and it depends only on the type of the voltage divider. Qgat fixes the 
amphfication performance of the PMT at a given high voltage, and it deter- 
mines the collectivity efficiency of the PMT. The amplification performance 
can be adjusted by the use of the appropriate voltage divider. 



3 Compcirisons of the model with the measurements 



3.1 The RI4O8 PMT case 



The 8-inch Hammamatsu R1408 PMT with 13 stages does not have a sin- 
gle photoelectron peak. Figure 5 shows the R1408 single photoelectron (spe) 
charge distribution. The main method to determine the gain G for this kind of 
PMTs is the photo-statistics method [5]. A series of LED/laser runs are taken 
at varying light levels, and for each run the statistical mean and variance of 
the spectra are measured as shown in Fig. 1, after that plot the variance versus 
the mean, and fit the linear region with a polynomial of the first order like in 
Fig. 2. The gain is extracted from the slope of the fit following the equation 
|51: 

^ = 2e«, (3) 
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Where E and o"^ are the mean and the variance of the ADC charge distribution, 
given by the following relations: 



E^q + lip (4) 
= 2qeG + (5) 

with e = 1.6 lO^^^C is the electron charge, q is the mean charge at the input 
of the ADC. Hp and are the mean and the variance of the ADC pedestal 
distribution. One important feature of the photo-statistics method is that it 
is ADC pedestal independent. 

The R1408 PMT was put under high voltage inside a dark box, illuminated 
by an LED driven by a fast pulser. The trigger output of the pulser was sent 
to Lecroy 222 Gate and Delay units, first to set a delay, and then to produce a 
200 ns gate which enclosed the PMT pulse that occurred when the LED was 
flashed, the using setup is shown in Fig. 3. The PMT signal was sent to the 
input of a Lecroy 2249W ADC, while the gate was coming from the 200 ns gate 
and delay output. The LED was flashed at 20 Hz at a variety of driving voltages 
from IV to 8V, all with very narrow (20 ns) driving pulse widths. In all cases, 
the PMT pulse was observed on an oscilloscope to be completely contained 
within the 200ns ADC gate. For each of the ADC distributions, the mean and 
variance were calculated, as it discussed above, and the ADC variances were 
plotted as a function of the ADC means which allowed determination of the 
gain Gq. The value of Gq was calculated from the hnear region obtained at 
low LED driving voltages. 

To compare the model to the data, the variance need to be written as function 
of the measured charge q. If the equation 5 is combined with the equation 1, 
we obtain the following relation: 

V laat 

Because the equation 6 is the theoretic variance, it does not contain any ADC 
pedestal variance term. 

For large q, the variance can be written like following: 

al^2eGoqltq'-^. (7) 

There are three cases of interest: 



(1) a < 2 : (7^ ^ oo. 
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(2) a = 2 : (7^ ^ 2eGog.at. 

(3) a > 2 : ^ 0. 

Figure 4 shows the variation of the variance obtained from equation 6 as func- 
tion of q, for different values of a. 



As shown in the Fig. 2 the variance shows a kind of saturation for high collected 
charge g, which corresponds to the case a = 2. 

Figure 6 shows the comparison of the model with the data obtained from 
the R1408 PMT at different applied high voltages, this PMT has a linear (or 
equally distributed) voltage divider. The reference [4] has more details about 
the voltage dividers and their different types, especially chapter 5. The value 
of Go depends on the applied high voltage, and it is determined as explained 
above from the slope of the hnear region. The value of Qsat is adjusted to 
reproduce the data, and the model reproduces the measurements with a — 2. 



3.2 Effect of the high voltage divider on the saturation level 



To improve the saturation performance of the PMT, a tapered voltage divider 
[4] is used instead of the linear one, but the gain is expected to become lower 
for the same total apphed high voltage. Details about the linear and tapered 
voltage divider, which were used during the measurements, will be given in a 
separate paper about the Double Chooz inner veto PMTs testing performed 
at Drexel University [5]. Figure 7 shows the comparison of the model with 
the data obtained from the measurements with the R1408 PMT at different 
applied high voltages, this time the PMT has a tapered voltage divider. As 
expected the tapered voltage divider increases the value of Qsat, which became 
at least 4 times larger. In addition. The model reproduces the measurements 
with a = 2. 



I want conclude this section with one remark about the effect of the voltage 
divider. Since the tapered divider enhances the collectivity efficiency, the gain 
can be expected to become lower at the same applied high voltage the gain 
due to the increasing of the collected charge at the anode. The assumption will 
be true if the gain depends only on the collected charge at the anode. However 
as it proposed by the equation 1, the gain does not depend on the collected 
charge q, but it depends on the ratio -2-. The tapered voltage divider will 

Qsat 

increase q and qsat at the same time, then if they are increased in the same 
way, their ratio will remain constant and the gain will not be different from 
the one obtained with the linear voltage divider at the same applied total high 
voltage. 
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3.3 The R5912 PMT case 



The 8-inch Hammamatsu R5912 PMT has a single photoelectron peak as 
shown in the Fig. 8. In order to apply the model to the case, a general for- 
mulation of the variance has to be found. The equation 5 can be written in a 
general way like following [5]: 

a' ^TgeG + crl. (8) 



7 is a factor depending on the single photoelectron charge distribution of the 
PMT. Since the interest of this paper is to explain the saturation profile of 
the PMT, it is enough to express the equation 6 hke following: 



where S = "jeGo, S is the slope of the linear region when the variance is plotted 
as function of the measured charge q. 

Figure 9 shows the comparison of the model with the data obtained from 
the R5912 PMT, this PMT came already encapsulated with its own voltage 
divider. The procedure for the gain measurement is exactly the same as for 
the R1408 PMT like described above. The model reproduces the data with 
q; = 2 in this case, too. 



34 The MCP-PMT case 



In this section, I compare the model to the data obtained with a new kind 
of PMT, the MCP-PMT. The ref. [4] contains a general description of this 
type of PMTs. The data which I used for the comparison with the model is 
published in the ref. [6]. To compare to the data, the following equation is 
used: 



G 



\ 



-1 + Wl + 8 



N, 



Ar2 



'ysat 



■ysat 



(10) 



The equation 10 gives the gain variation as function of the number of the 
incident photons N^. N^sat is the number of incident photons at which the 
PMT saturation occurs. N,ysat is directly determined from the data as the 
value of which corresponds to the ratio ^ = 75- The demonstration of 
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the equation 10 is made in the appendix. Figure 10 shows the results of the 
comparison. Also in this case, the model reproduces the data. 



4 CONCLUSION 

To describe some of the known properties of the photomultiplier tubes, a PMT 
model is proposed. The PMT characteristics like the quantum efficiency, the 
collectivity efficiency and the gain are linked together in natural way as demon- 
strated in the appendix. The model is successful to describe the saturation 
profile measured for three different PMTs with different single photoelectron 
charge distributions. The comparisons with the measurements show that the 
value of the coefficient a is 2. The effect of the voltage divider - on the gain 
variation and the saturation performance - is explained in the model by the 
variation of the ratio — . 

Qsat 
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APPENDIX: Formulation of the PMT characteristics 

The collected Charge at the anode and the gain: 

The collected charge at the anode can be expressed as function of the number 
of the photoelectrons, which are collected at the first dynode and the gain 
G as following: 



Where e = 1.6 lO^^^C is the electron charge. 

Now, if the equation 1 with a = 2 is taken in account, the following equation 
will be obtained: 



q = iJ,Ge, 



(A-1) 




G^ + G^-Gl^ 0. 



(A-2) 
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The equation A-2 can be considered as the master equation of the PMT, which 
hnks the observables of the PMT to each other, the quantum efficiency via 
the collectivity efficiency represented by Qgat and the gain G. 

The equation A-2 has hke solution the following expression: 



eG 



\ 



'isat 



laat 



(A-3) 



The equations A-1 and A-3 let to express the collected charge at the anode, 
like following: 



2 

Isat 



'isat 



(A-4) 



In the case of the cathode current saturation, jj, becomes constant which makes 
the collected charge at the anode and the gain constants, too. 

The quantum efficiency: 

If the equations 1 and A-1 are put together, /j, can be expressed as function of 
the collected charge at the anode, q, as following: 




(A-5) 



IJ, is linked to the number of incident photons, Nj, by the following relation: 
1^ = p7]N^, (A-6) 



Where p is the collectivity efficiency between the photocathode and the first 
dynode, r] is the quantum efficiency of the PMT. 

The equations A-5 and A-6 let to express the product prj as following: 
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The number of photoelectrons at the saturation: 



For q = Qsat, the equation A-5 gives the number of photoelectrons at the 
saturation, figat, collected at the first dynode: 

f^sat = V2^, (A-8) 
eLrn 



Usat characterizes the saturation level of the PMT at a given high voltage. If 
the equation 2 is used with the equation A-8, fisat can be expressed as function 
of the applied high voltage: 



Since the value of Qsat depends only on the type of the voltage divider, the 
number of the photoelectrons at the saturation decreases with increasing ap- 
plied high voltage for a given PMT with a given voltage divider. 
Using the equations A-3, A-6 and A-8, the gain can be written as function of 
the number of incident photons A^-^: 



\ 



-1 + .1 + 



A "t 

^sat 



(A-10) 



Where N^^at = —■ 
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R1408 PMT ADC charge distribution 




ADC count 



Fig. 1. ADC distribution of the PMT pulses for LED flashing with 3V from pulser. 
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Fig. 2. ADC variance((7^) versus ADC mean(g) for a variety pulse heights. The 
slope of the fit of the linear region gives a gain/channel=slope/2=24.9 chan- 
nel / photoelectron. 
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Fig. 3. The test setup scheme. 
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Fig. 4. The ratio -^j^ as function of the mean of the cohected charge at the anode 
q, as given by the equation 6 with qgat = 500. 
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R1408 single photoelectron charge distribution 




Fig. 5. The single photoelectron (spe) charge distribution for the 8-inch, 13 stages, 
R1408 PMT. The ADC pedestal is not removed. The signal from the PMT output 
is amplified 10 times. 
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R1408 PMT with linear voltage divider 
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Fig. 6. Comparison of the model with the data obtained at different high voltages, 
Qsat and eGo are in ADC channels. The stars represent the data, and the dashed 
line represents the model with a = 2 for all. 
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R1408 PMT with tapered voltage divider 
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R1408 PMT with tapered voltage divider 
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Fig. 7. Comparison of the model with the data obtained at different high voltages, 
Qsat and eGo are in ADC channels. The stars represent the data, and the dashed 
line represents the model with a = 2 for all. 
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R5912 single photoelectron charge distribution 
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Fig. 8. The single photoelectron (spe) charge distribution for the 8-inch R5912 PMT. 
The ADC pedestal is not removed. The signal from the PMT output is amplified 
10 times. 
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Fig. 9. Comparison of the model with the data obtained at different high voltages, 
Qsat and S are in ADC channels. The stars represent the data, and the dashed line 
represents the model with a = 2 for all. 
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MCP-PMT 
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Fig. 10. Relative gain ^ vs photon rate N^. The dashed hne represents the model 
(i.e. the equation 10 with N^sat = 2.5 lO^count/cm^/s), and the stars represent the 
data from the ref. [6] . 
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